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1. Introduction 

B Physics, and in particular the study of B and B s meson decays and mixing, remains an 
important part of Flavor Physics. Such studies enable consistency tests of the Standard Model and 
placing of bounds on New Physics effects. Lattice QCD is playing a crucial role in this effort by 
providing the necessary nonperturbative QCD inputs. 

This talk describes recent results by the HPQCD collaboration on B meson mixing parameters, 

/b, \/Bb s , fB d \]BB d and their ratio t, , the first fully consistent Nj ■ = 2 + 1 lattice QCD calculations 
of these quantities [[I]]. We work with four of the MILC coarse (a « 0.12fm) lattices and two of the 
fine (a m 0.09fm) lattices. We use NRQCD for the Z?-quark and the AsqTad action for both valence 
and sea light quarks. Details are given in [[I]]. In the Standard Model the mass difference in the 
Bq — Bq (q = d,s) system is given by 



AM, 



w t 



\V tq v; b \ 2 ri%S (x t )M Bg f£B B . 



2 D 



6tt 2 



(1.1) 



In addition to well determined quantities, eq.( |l.l[ ) involves the combinations of CKM matrix ele- 
ments \V tq V* h \ 2 and the nonperturbative QCD factors fg B Bq . The latter are determined from the 
matrix element of the four-quark operator are color indices that are summed over), 



OL = 



%(V-A)% Vi(V-A)yJ q 



sandwiched between the B q and B q states. 



(OL) (n) = <5,|0L| W) ee -ti q B Bq (/i)Mi 



(1.2) 



(1.3) 



One sees from eq.(|l.l[) that the ratio of the two CKM matrix elements \V t d\ and |V (S | can be deter- 
mined from 

\V td \ e lAM d M Ba 



\V a 



once theory provides the important ratio 



AM, M Brl 



fB s \/B~B~ s 
fB d \fB~B~ d 



(1.4) 



(1.5) 



One of the main goals of [|l|] was to produce a state-of-the-art lattice calculation of E, . We were able 
to determine this quantity with total error of 2.6%. 



2. Hadronic Matrix Elements in the Effective Theory 

The b quarks in our calculations are described by an effective theory NRQCD. Instead of 



the full QCD field of eq.(jl.2|) one works with fields *Pg that create a heavy quark or *Pg that 
annihilate a heavy anti-quark. The [V — A] x [V — A] four-fermion operator becomes, 



OL eff EE 



%(V-A)% *fyy-A)% + Vq(V-A)% %(V-A)V q 



\TJJ 



(2.1) 
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In the effective theory one finds mixing with another operator at one-loop of [S — P] 
structure, even at lowest order in 1 /M. 



OS eff = 



[S-P] 
(2.2) 



Furthermore there is a tree-level dimension 7 correction to OL e ^ at G{\jM). 
1 r/ W e -7(V-A)%)(%(V-A)^ 



OLjl = 

J 2M L 



+ 



(2.3) 



+ (W Q (V -A)W q ) [VWq ■ y{V-A)^ q ^ 

We have calculated the matching between matrix elements of the effective theory operators Olff f, 
OS ef f and OLjl and (OL) of eq.Q in full QCD through &(a K ,A QCD /M, a s /{aM)) and find 

(OL) m (ii) = [\ + a s pn}(OL e ff) + a sPl2 (OS e ff} + 

(OLjl) - a, [ C 11 (OlfH) + C 12 (OS'H) ] + 0{al olAqcd/M). (2.4) 

The matching coefficients pn, p\ 2 , £ n an d £ 12 we listed (for /i =Mb) in [||]. 

The goal of lattice simulations is to obtain the matrix elements (O), with O = Olfff, OS e ^ or 
OLjl. To this end one calculates three-point correlators, 



C%\t u t 2 ) = £ <0|*g (xuh) [a 6 O(0)} ^ 2 ,-t 2 )\0) 



X U %2 



together with two-point correlators 

-Zpt 



c 2 a p)= I <o|*2,(*,0 <(^2,o)|o>. 



(2.5) 



(2.6) 



X\ ,X2 



One works with dimensionless operators a 6 which are kept fixed at the origin of the lattice. <J>^ 
is an interpolating operator for the B q meson of smearing type "a", and spatial sums over x\ and 
x 2 ensure one is dealing with zero momentum B q and B q incoming and outgoing states. Eq.(2.5) 
corresponds to a B q meson (or excitations thereof with the same quantum numbers) being created 
at time — t 2 which then propagates to time where it mixes onto its anti particle a B q . The latter 
meson is annihilated at time t\ . 

We have accumulated simulation data for a range in t\ , t 2 for three-point and in t for two-point 
correlators, i.e. 1 <h,t 2 ,t < T max with T max = 24 on the coarse lattices and T max = 32 on the fine 
ensembles. As is well known, staggered quarks lead to meson two-point correlators that have both 
regular and time oscillating contributions. Hence a fit ansatz for eq^ZjS) would be 

N-l 



N-l 

t V uatP „-4(f-l) 



7=0 k=0 



-i)T 



(2.7) 



Similarly, the presence of oscillatory contributions makes fitting the three-point correlators, eq.([2.5|) 
particularly challenging. The appropriate ansatz is 



N-lN-1 

E L A 

7=0 k=0 



JV— 1JV— 1 

-Ej(ti-1) e -E k {h-\) + y y B aP ^_jy, g -B;(fi-l) e -E k {t 2 -\) 
7=0 k=0 



i OA (I 



N-lN-l 
7=0 k=0 



(2.8) 
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In terms of the fit parameters in (2.7) and (2.8) the matrix elements we are interested in and which 



appear on the RHS of ( |2.4| ) are given by, 



_ K 2M B An,f 

(B q \OB q )= *« °° (2.9) 



In order to extract -^j, we have carried out simultaneous fits to C^P(t\,t 2 ) f° r ce = j8 together 

with a matrix of two-point correlators CjZ (?) with all possible combinations of smearings at source 
and sink. We used Bayesian fitting methods and typically employed N = N = 4 ~ 6 number of 
exponentials. 

3. Chiral and Continuum Extrapolations 

After the fits to simulation data described above (see [|J] for more details) one ends up with 
lattice results for rV 2 f B „\ M B B B for each of the 6 MILC ensembles that we worked with. Here 
r\ is a scale derived from the static potential which can be used to make dimensionful quantities 
dimensionless. We exploit the fact that the MILC collaboration has calculated r\ja for each of 
their ensembles. In order to make contact with the real world one needs to extrapolate both to the 
continuum (lattice spacing a — > 0) and to the chiral (mu g h t — ► m u u) limits. We do so by fitting to 
the following ansatz, inspired by chiral perturbation theory (ChPT). 



fB q \JM B B Bq =ci [l + -A/ g + c 2 (2m f + m s )n +c 3 m (? r 1 ] x [1 +c A a s {a/r l ) 1 + c 5 (a/n) 4 ]. 

(3.1) 

Af q includes the chiral logarithms including those specific to Staggered ChPT and was calcu- 
lated by C.Bernard, J.Laiho and R.Van de Water We show chiral/continuum extrapolations 
for r\ ,2 f Bs ^M B B Bs and rf Z f Bd ^M Bi B Bi in Fig.l and for % yj M Bs / M Bd and for B B$ in Fig.2. 
The red curves are the continuum exttapolated curves and the red circle gives the results at the 
physical point. Table I gives our error budget. 

4. Results 

Using central values coming from the physical (red) points in the figures and the errors summarized 
in Table I, we can now present our main results [jl|]. 

*- h ^ Bl) l.258(25)(21), B Bs = 1.33(5) (3) (4.1) 



fB d \JB Bd 
and using r\ = 0.321 (5)fm [||], 

0.321 \ 3/2 „ p— _ N /0.32O 3/2 



/ Bj V^=266(6)(17) MeV, f BdX J B Bd = 216(9) (12) j MeV, 

(4.2) 

where the first error comes from statistics + chiral extrapolation and the second is the sum of all 



other systematic errors added in quadrature. The result for f B \ B Bs in eq.(4.2) is consistent with 
but more accurate than our previously published value of 281(21)MeV [¥J|. 
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Figure 1: Chiral and continuum extrapolation of jy /b,\/ Mb s Bb s , r x ' fB d \ Ms d BB d . The red curve shows 
the continuum extrapolation and the "physical point is at mn g i, r /m s = 1 /27.4. 



Combining our lattice result for E, with the experimentally measured mass differences AM^ 
0.507 ±0.005 ps~ l @ and AM, = 17.77 ±0.10 ±0.07 ps' 1 [§] leads to, 



m 

\Vts\ 



0.214(1)(5) 



(4.3) 



where the first error is experimental and the second from the lattice calculation presented here. Our 
Bag parameter result Bg ! can be combined with experimental AM S and t(B s ) to form the Standard 
Model prediction for the branching fraction Br(B s —* /i + ;U~) [PI]. We find 



Br(B s -» pL + ll~) = 3.19(19) x 10" 
improving on the previous accuracy available. 



(4.4) 



source of error 




fB d \J B Bd 


5 


stat. + chiral extrap. 


2.3 


4.1 


2.0 


residual a 2 extrap. 


3.0 


2.0 


0.3 


uncertainty 








r\l 2 uncertainty 


2.3 


2.3 




8b*bk uncertainty 


1.0 


1.0 


1.0 


m s and m/j tuning 


1.5 


1.0 


1.0 


operator matching 


4.0 


4.0 


0.7 


relativistic corr. 


2.5 


2.5 


0.4 


Total 


6.7 


7.1 


2.6 



Table 1: Errors in % for f Bs J B Bs , fB d y BB d and | . 
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Figure 2: Same as Fig.l for % yjM B jM Bi and for B Bs 



Our fits to a matrix of B q two-point correlators has also allowed us to update previous results on B c i 
and B s decay constants. The final numbers including all errors added in quadrature become, 



^ = 1.226(26), 



f Bd = 190(13) 



0.321 

n \fm] 



3/2 



MeV, f Bs = 231(15) 



0.321 

n [fm] 



3/2 



MeV. 



(4.5) 



(4.6) 



These results for f Bq are consistent with but about one a lower than the values f Bd = 216(22)MeV 
and f Bs = 260(29)MeV given in [ 1C, 11]. The main difference between the analysis carried out here 
and in Jl0| ] is that in the latter case chiral extrapolations were done based only on coarse lattice data 
and furthermore no attempt was made to extrapolate explicitly to the continuum limit. 



5. B s Meson Decay Constant Using Other Actions 

The HPQCD collaboration has initiated a project of studying B physics with NRQCD ft-quarks 



and HISQ light quarks [12]. The hope is that by going from the AsqTad to the more highly im- 



proved HISQ action discretization errors will be further reduced, in particular those coming from 
taste breaking effects. We are also simulating the B c system with NRQCD b- and HISQ charm- 
quarks [|l3|]. Here we present first preliminary results for the B s meson decay constant f Bs in Fig. 3. 
One sees that the slope versus lattice spacing is significantly reduced as one moves to the HISQ 
action. Furthermore both actions lead to the same continuum limit within errors. We are pursuing 
other approaches to B physics as well, in particular one based on relativistic HISQ fe-quarks. Initial 
calculations using relativistic ^-quarks show results for f Bs in good agreement with the NRQCD- 



AsqTad and NRQCD-HISQ values presented here [|1J. 
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Figure 3: Comparison between NRQCD-AsqTad and NRQCD-HISQ calculations of the B s meson decay 
constant. 
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